Elastic electron scattering in quantum corrals: The importance of the shape of the 

adatom potential. 
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We report elastic scattering theory for surface electron waves in quantum corrals defined by 
adatoms on the surface of noble metals. We develop a scattering-matrix technique that allows us to 
account for a realistic smooth potential profile of the scattering centers. Our calculations reproduce 
quantitatively all the experimental observations, which is in contrast to previous theories (treating 
the adatoms as point scatterers) that require additional inelastic channels of scattering into the 
bulk in order to achieve the agreement with the experiment. Our findings thus indicate that surface 
states are not coupled to the bulk electrons. 
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Advances in modern nanotechnology made it possible 
to manipulate adatoms on the surface of a metal, arrang- 
ing them into ordered structures coined as "quantum cor- 
rals". Using the scanning tunneling microscopy (STM), 
Crommie et aZi Vi 3 studied the scattering of the surface 
electron waves residing at (111) faces of Cu. These sur- 
face states interact strongly with Fe adatoms, and the 
spatial variation of the STM differential conductance re- 
vealed beautiful images of the standing wave patterns in 
the quantum corrals. In addition, the experiment showed 
a series of pronounced peaks in the energy spectrum of 
the differential conductance dl/dV at the center of the 
corrals. 

In order to describe the experimental observation^, 
Heller et alA have developed the multiple-scattering the- 
ory for surface electron waves in quantum corrals. In 
their theory each adatom was modelled as a point-like 5- 
function potential supporting isotropic s-wave scattering. 
The quantitative agreement with the experiment was 
achieved by assuming an additional (inelastic) channel 
of scattering (presumably into the bulk metallic states) . 
The authors concluded that absorption is the dominant 
mechanism for the broadening of the energy levels seen 
in the experiment, and estimated that ~ 25% of the in- 
cident amplitude is reflected, ~ 25% is transmitted, and 
~ 50% is absorbed. The importance of the electron scat- 
tering to the bulk states for the level width broadening 
in the quantum corrals was also asserted by Crampin et 
al^ and Cramplin and Bryant. 

An alternative purely elastic scattering theory for the 
same quantum corral structures was reported by Harbury 
and PorodZ. They modelled the adatoms by finite-height 
potential barriers, as opposed to "black dot" 5-point ab- 
sorbing scattering potential adopted in the above cited 
works^^. The elastic theory accounts well for the spa- 
tial variation of the wave function in the quantum corrals, 
but overestimates the broadening of the resonant levels, 
especially for higher energies. The findings of Harbury 
and Porod therefore suggest that the features of the spec- 
trum can be sensitive to the detailed shape of the scat- 



tering potential. 

It is important to stress that accounting for a de- 
tailed shape of a scattering potential was crucial for 
quantitative description of many phenomena in quantum 
nanostructures. This for example includes the Hall and 
bend resistance anomalies in four-terminal junctions 8 , 
the breakdown of quantized conductance in quantum 
point contacts calculated using realistic potentials^, the 
explanation of a branched flow in a two-dimensional elec- 
tron gas 10 . In the present Brief report we develop a scat- 
tering matrix approach that allows us to account for a re- 
alistic smooth potential of the adatoms. We demonstrate 
that for such a potential the broadening and positions of 
the resonant states as well as the scattering wave function 
in the quantum corrals can be quantitatively described 
by the inelastic theory alone, without the assumption of 
any additional (inelastic) scattering channels. Our find- 
ings thus support the conclusions of Harbury and Porodi 
that the adatoms do not significantly couple the surface- 
state electrons to the bulk metallic states. 

The differential conductance dl/dV of the STM tun- 
nel junction is proportional to the local density of states 
(LDOS) which is given in terms of the scattering eigen- 
states of the Hamiltonian H, ijj q (r)&, 

dl/dV ~ LDOS(r, E) = £ \yb q {v)\ 2 5{E - E q ). (1) 



In order to calculate the scattering eigenstates ip q (r) 
we adopt to the problem at hand the scattering matrix 
tcchniqueii that was recently developed for numerical 
solution of the Helmholtz equation for resonant states 
of dielectric optical cavities with both complex geometry 
and variable refraction index. This is possible because 
of a direct correspondence between the Helmholtz and 
Schrodinger equations 12 . The advantage of the scatter- 
ing matrix technique is that it provides an efficient way 
to treat the smooth realistic profile of the adatom. Note 
that commonly used methods based on the discretization 
of the scattering domain would be rather impractical in 
terms of both computation power and memory, because 
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the smoothly varying potential of the adatom has to be 
mapped into a discrete grid with a very small lattice con- 
stant. 

We consider a two-dimensional ring-shaped corral 
structure. Experimental observations suggest that 
adatoms strongly disturb the local charge density at the 
finite distance ~ 7Ai An exact experimental shape of 
the adatom potential is not available. We thus model 
this potential as a Gaussian with the half-width a and 
the height V centered at the location (x ,y o ), V(x,y) — 
V exp[— (x — x ) 2 /2a 2 ] exp[— (y — y ) 2 /2a 2 ], (see below, 
inset to Fig. EJ. 

In order to calculate the scattering eigenstates we di- 
vide the quantum corral in inner, outer and intermediate 
regions. In the inner region (inside the corral) and in the 
outer region (outside the corral) the adatom potential is 
negligible, V(x, y) = 0. Therefore, in these two regions 
the solution to the Schrodinger equation can be written 
in analytical form. Introducing the polar coordinates we 
can write for the wave function outside the corral 

+00 

yout = [A q H^{kr)+B q H^{kr))e^, (2) 

q— — oc 

where H q 2 ^ are the Hankel functions of the first 

and second kind of the order q describing respectively 
incoming and outgoing waves; k — \/2m*E/fr, with m* 
being the effective electron mass. The expression for the 
wave function inside the corral \& m can be written in a 
similar fashion as an expansion over Bessel functions J q . 

We introduce the scattering matrix S in a standard 
fashioniiiiS, B — SA, where A, B are column vectors 
composed of the expansion coefficients A q , B q for incom- 
ing and outgoing states in Eq. @. The matrix element 
S q iq gives the probability amplitude of scattering from 
an incoming state q into an outgoing state q' . In order to 
apply the scattering matrix technique we divide the inter- 
mediate region (i.e. the region where the adatom poten- 
tial V(x,y) is distinct from zero) into N narrow concen- 
tric rings. At each i-th boundary between the rings we 
introduce the scattering matrix S 1 that relates the states 
propagating (or decaying) towards the boundary, with 
those propagating (or decaying) away from the bound- 
ary. The matrices S 1 are derived using the requirement 
of continuity of the wave function and its first deriva- 
tive at the boundary between the two neighboring rings. 
Successively combining the scattering matrixes for all the 
boundarie s- 1 -™' ~ 2 , S 1 (£> . . . (£> S N , we can relate the com- 
bined matrix to the scattering matrix S. With the help 
of the scattering matrix S we determine the wave func- 
tion in the outer region for every incoming state q in Eq. 
@. Using the expression for the matrixes S 1 we then 
recover the corresponding wave functions in the interme- 
diate region as well as the wave function ^> m in the inner 
region. 

Note that in the scattering matrix technique one com- 
bines only two scattering matrixes at each step. Hence, 
it is not necessary to keep track of the solution for the 



wave function in the whole space. This obviously elim- 
inates the need for storing large matrices and facilitates 
the computational speed. 

Using our scattering matrix technique we calculate the 
bias voltage dependence and the spatial distribution of 
the LDOS for 60-Fe-adatom, 88.7-A-radius circular quan- 
tum corrals reported by Heller et a/4(Figs. Q]and|2J). The 
Fe-adatoms are placed on the meshes of a 2.55A trian- 
gular grid corresponding to the hexagonal Cu(lll) lat- 
tice. The effective mass used in all the simulations was 
taken m* = 0.361mo and the electron band-edge energy 
Eq = 0.43eV below the Fermi energy of the electronsii 
In the absence of applied voltage V these parameters cor- 
respond to the wavelength of electrons A = 30A . For the 
parameters of the adatom potential we usa£ Vq = 2.5eV, 
a = 1.52A which correspond to those used by Harbury 
and Porodi who modelled the adatoms as hard wall finite 
potential barriers of 1.52-A diameter of the height of 2.5 
eV. 

Figures also shows corresponding results of the mul- 
tiple scattering theory of Heller et al. 4 . Both theories 
show a similar level of agreement with the experimental 
data in both the peak positions and level broadcnings 
for the differential conductance as well as in the number 
and the peak positions for the spatial LDOS distribution 
throughout the cross-section of the quantum corral. But 
in contrast to the case of 5-scatterers used in Ref. our 
model agrees quantitatively well with the experimental 
data without introduction additional inelastic scattering 
channels. 
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FIG. 1: The experimental spectrum of the differential con- 
ductance dl/dV at the center of the 88.7-A-radius 60-Fe- 
adatom circular quantum corral structure on Cu(lll) sub- 
strate (adopted from Ref. 0) (solid curve). The calcu- 
lated spectrum, dotted line: Our scattering matrix technique 
applied for a smooth adatom potential with Vb = 2.5eV, 
a = 1.52A; dashed line: Multiple-scattering theory for S- 
barrier adatom potential with inelastic channel of scattering 
(adopted from Ref. 0). 

Figure |3 represents the differential conductance spec- 
tra dl/dV at the center of the quantum corral struc- 
ture for various widths of adatom potentials a. The in- 
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FIG. 2: The experimental curves (solid lines, adapted from 
Ref. 0) for the LDOS subject to the tip position inside a 
circular corral for low bias voltages. The calculated LDOS, 
dashed line: Our scattering matrix technique applied for a 
smooth adatom potential with Vb = 2.5eV , a = 1.52A; dot- 
dashed line: Multiple-scattering theory for 5-barrier adatom 
potential with inelastic channel of scattering (adopted from 
Ref. 0). Parameters of the structure and are the same as 
those in Figs. Voltages are given in volts and measured 
relatively to the bottom of the surface-state band. All the 
theory voltages are shifted by -0.01 V relatively to the exper- 
iment. 



crease in the potential width leads to significant narrow- 
ing of the peaks in the spectra due to stronger confine- 
ment of the wave function inside the corral. The best 
agreement with the experiment is achieved for a ~ 1.5A, 
which is in agrement with the results reported by Har- 
bury and Porodi It is often assumed that because the 
spatial extent of the scattering potential (a ~ 7 A) is 
small compared to the wavelength of the incoming par- 
ticles (A ~ 30 A), the adatom potential can be treated 
as a point scatterer or even as a continuous boundary 
V(r) ~ S(r — ro)^*i^. The results presented in Fig. |3| 
clearly show that even though a <C A, the finite width 
of the scattering potential affects strongly the observed 
characteristic of the systems. Our calculations thus sig- 
nify the importance of the shape of the scattering po- 
tential for achieving the quantitative agreement with the 
experiment of the voltage dependence of dl/dV. 

The experiment suggests that for high energies of in- 
coming electrons (i.e. for large tip voltages) the signifi- 
cant fraction of Fe ad atoms can move from their original 
positions. We therefore study the effect of this displace- 
ment on the shape and broadening of the resonant states 
of the quantum corral. Figure0|shows the calculated dif- 
ferential conductance for a quantum corral for the case 
when the adatoms of an ideal circular corral are randomly 
shifted from their positions by one and two lattice con- 
stant (as illustrated in the insets to Fig. HJ. One-lattice 
shift does not seem to have a significant effect on the 
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FIG. 3: The calculated spectrum of the differential conduc- 
tance dl/dV at the center of the circular quantum corral 
structure for various widths a of adatom potential for fixed 
barrier height Vo = 2.5eV. Parameters of the structure are 
the same as those in Fig. Q Inset shows the schematic shape 
of the potential. 
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FIG. 4: Effect of the non-ideal positioning of the adatoms on 
the differential conductance of the quantum corral. Parame- 
ters of the structure and the adatom potential are the same 
as those in Figs. 0and|5| Inset illustrates the displacement 
of the scatterers from their ideal positions (white circles) in a 
circular geometry for two-lattice shift. 



broadening of the peaks in the differential conductance. 
The deviation from an ideal circular case become rather 
noticeable for the shift in adatom positions by two lattice 
sites. As expected, these deviations are more pronounced 
for larger energies of incoming electrons. Note that we 
performed simulations for different realization of the en- 
sembles of scatterers (keeping the average shift fixed to 
one or two lattice sites) , and all of them are almost indis- 
tinguishable. This is because of the self-averaging char- 
acter of scattering in a quantum corral due to a large 
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number of scatterers. Our calculations thus pinpoint the 
deviation from a regular arrangements of scatterers as 
an additional factor that should be taken into account 
for achieving the quantitative agreement with the exper- 
iment. 

To conclude, we studied the scattering of electron 
waves by quantum corral structure for the case of realis- 
tic smooth potential of adatom scatterers. We achieved 
a detailed agreement with the experiment without in- 
troducing the additional inelastic channels of scattering. 



Our results thus support the theory that adatoms do not 
significantly couple the surface-state electrons to the bulk 
states. Our findings also suggest that accounting for a re- 
alistic potential shape may be of particular importance 
for the quantitative description of quantum mirages ex- 
periments. 
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